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Very short quantum wires (quantum contacts) exhibit a conductance structure at a value of 
conductance close to 0.7 x 2e 2 /h. It is believed that the structure arises due to the electron- 
electron interaction, and it is also related to electron spin. However details of the mechanism of 
the structure are not quite clear. Previously we approached the problem within the restricted 
Hartree-Fock approximation. This calculation demonstrated a structure similar to that observed 
experimentally. In the present work we perform restricted and unrestricted Hartree-Fock calcula- 
tions to analyze the validity of the approximations. We also consider dependence of the effect on 
the electron density in leads. The unrestricted Hartree-Fock method allows us to analyze trapping 
of the single electron within the contact. Such trapping would result in the Kondo model for the 
"0.7 structure". The present calculation confirms the spin-dependent bound state picture and 
does not confirm the Kondo model scenario. 

PACS: 73.61.-r, 73.23.Ad, 71.45.Lr 

The quantized conductance G — nG2, n — 1, 2,3, G*2 = 2e 2 /h 1 through a narrow quantum point contact was 
discovered in 1988 1 ' 2 . This quantization can be understood within a one-dimensional (ID) non-interacting electron 
gas picture, see e.g. Ref. 3 . In the present work we are interested in a deviation from the integer quantization. 
This deviation, the so called "0.7 structure" has been found in experimental works ' 5 . The structure is a shoulder- 
like feature or a narrow plateau at G ss 0.7G2- More recent work demonstrates that there are some above barrier 
excitations related to the structure 6 , and that the structure evolves down to G w O.5G2 in longer quantum contacts 7 . 
Dependence of the structure on the longitudinal magnetic field has been studied already in the pioneering work 4 . 
This study clearly demonstrated that the effect is somehow related to the electron spin. Authors of a recent 
experimental work 8 argue that the structure signals formation of a Kondo-like correlated spin state. 

There have been suggestions to explain the "0.7 structure" by spontaneous magnetization of the ID quantum 
wire 9-13 , or by formation of a two-electron bound state with nonzero total spin 14,15 . These suggestions implicitly 
assume that 2D leads connected to the contact are qualitatively important for the effect because there is the 
rigorous Lieb-Mattis theorem 16 that claims that the ground state of a ID many-body system has zero spin. 

A Hatree-Fock calculation of the conductance has been performed in the Ref. 17 . This calculation demonstrated 
a structure similar to that observed experimentally. The ground state has zero spin in accordance with the Lieb- 
Mattis theorem, but nevertheless the structure found in the calculation is intrinsically related to the spin because 
it disappears without account of the exchange electron-electron Coulomb interaction. The structure is related to 
the formation of the charge density wave within the contact or in other words to the spin-dependent bound state 
within the contact. 

The present work has been stimulated by the recent suggestion that the 0.7-structure signals formation of a 
Kondo-like correlated spin state 8 , see also Ref. 18 . The restricted Hartree-Fock (RHF) approximation employed in 
Ref. 17 is not sufficient to follow this suggestion. However the unrestricted Hartree-Fock (UHF) approach can shed 
light on the problem. In the present work we consider only zero temperature case. The RHF method implies that 
spin up and spin down single electron orbitals are identical while in the UHF method those orbitals are completely 
independent. The RHF is explicitly rotationally invariant, but it is not very effective in accounting for electron- 
electron correlations. The UHF is much better at accounting for the correlations, but it violates the rotational 
invariance. There is no doubt that even UHF cannot account for the long range Kondo-like dynamics. However 
it can indicate localization of a single electron within the contact. This would immediately imply the Kondo-like 
dynamics. Our calculation shows that such localization can take place in longer contacts and at low electron 
density in leads. However it always leads to a very special dependence of conductance on the gate voltage which 
is different from that observed experimentally. In the regime when the dependence of conductance on the gate 
voltage is similar to the experimental one the results of RHF and UHF are practically identical and this indicates 
the validity of both approximations. We also study dependence of the "0.7 structure" on the electron density in 
the leads. The structure disappears at high density and it is getting more pronounced at the low density in a 
qualitative agreement with experiment. 

It is well known, see Ref. 3 , that in the independent particle approximation, i.e. in the case of an ideal electron 
gas, the conductance for a given transverse channel is proportional to the barrier transmission probability at Fermi 
energy T, 

2e 2 

G = fL.T. (1) 
In case of interacting particles this formula should be also valid because before and after the potential barrier 



the density of electrons is high enough, and hence the interaction is negligible. However one cannot use a single 
particle description to calculate the transmission probability T because in the vicinity of the barrier the electron 
density is low, and hence the many-body effects are very important. To calculate the transmission probability T 
the following method is applied. Consider electrons on a ID ring of the length L with a potential barrier of the 
length I somewhere on the ring. It is important that L 3> I. There is no current in the ground state of the system. 
Now let us apply a magnetic flux through the ring. This flux induces the electric current. Note that it is not a 
real magnetic field, this is a fictitious gauge field that generates the current without applying any voltage. It is 
especially convenient to take the gauge field that provides the Bohm-Aharonov phase ip = ir/2. We use this choice 
in our calculations. The induced current can be calculated by solving many-body Schrocdingcr equation. It can 
be an exact solution or an approximate one like RHF or UHF. It has been demonstrated in Ref. 17 that to find the 
barrier transmission probability at Fermi energy one has to solve the many-body problem twice: without the barrier 
and with the barrier. The ratio of electric currents squared gives the transmission probability, T — (Ju/Jo) 2 - This 
formula is valid without an external magnetic field. Repeating considerations of Ref. 17 one can prove that with 
the magnetic field, i.e. with the spin splitting, the effective transmission probability is given by 
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where is the electric current of electrons with spin up and spin down correspondingly. Equation (2) can 
be also applied for UHF calculations without an external magnetic field. The relation T = (Ju/Jo) 2 has been 
applied recently to study conductance through a system of strongly correlated spinless fermions 19 . In this work 
the many-body problem has been treated exactly via the Density Matrix Renormalization Group algorithm. 
The Hamiltonian of the many body system we consider is of the form 
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where indexes i and j numerate electrons, Xi is the periodic coordinate on the ring of length L (0 < x < L), 
A = n/2L is the fictitious gauge field. The electron-electron Coulomb repulsion is of the form 
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where a t « 2 is the effective width of the transverse channel, see Ref. 17 , and D(x, y) is the length of the shortest arc 
between the points x and y on the ring. We use atomic units, so distances are measured in unites of Bohr radius, 
as = eh 2 /me 2 , and energies are measured in units of E unit = me 4 /Ji 2 e 2 , where m is the effective electron mass and 
e is the dielectric constant. For experimental conditions of works 4-8 these values are the following: ~ 10~ 2 /ito, 
Eunit ~ 10~ 2 eV. To model the gate potential we use the following formula for the potential barrier 



U(x) = 

Plots of U(x) for 1=8,10,12 are shown in Fig.l. 
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FIG. 1. The "gate" potential (5) at U = 4 and I = 8, 10, 12. 



To solve the many-body problem described by the Hamiltonian (3) we use the Hartree-Fock (HF) approximation. 
In the HF approximation the many body wave function is represented in the form of the Slater determinant of 
single particle orbitals Lp ia {x). The index i shows the coordinate state of the orbital, and the index a = ±1/2 shows 
the spin state of the orbital. Each orbital obeys the equation 



h<pi. 



(6) 



where ti a is the single particle energy and h is the HF Hamiltonian 
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The summations are performed over all filled orbitals. In the Restricted Hartree-Fock (RHF) method an additional 
constraint, ip^{x) = ifi^(x), is imposed on the solutions of Eqs. (7). This provides rotational invariance of the 
solution. In the Unrestricted Hartree-Fock method (UHF) the additional constraint is omitted. As a result the 
UHF method provides much better account of electronic correlations. The price for this is a spontaneous violation 
of the rotational invariance. 

For computations we use a finite grid. In the grid modification of the Hamiltonian (7) the kinetic energy (p— A) 2 ip 
is replaced by [2|cp(n) | 2 - ip*(n + l)e lAh tp(n) - tp*(n)e~ lAh tp(n + 1)] /2h 2 . Here h is the spacing of the grid and 
ip(n) is the wave function on site n of the grid. The electric current corresponding to the grid Hamiltonian reads 
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The current is conserved because of the gauge invariance of HF equations. 

For computations we use a grid of 400 points on a ring of length L = 80. Total z-projection of the spin is zero, 
so the number of electrons with spin up is equal to that with spin down, N]' — A^. We perform calculations for 
the total number of electrons N = + Aj_ = 78, 118, 158. This corresponds to the following values of the number 
density of electrons on the ring: n = N/L w 1, 1.5,2. This is the effective linear density, therefore one cannot 
compare n quantitatively with the density of electrons in real two-dimensional leads used in experiments 4-8 . 
However a qualitative comparison is possible: the smaller the real density, the smaller no, and hence the smaller 
the Coulomb screening. Results of calculations for three different values of the barrier length, I = 8, I = 10, I = 12, 
and three different values of the electron density in the "leads", no ~ 1, n o ~ 1-5, and n ~ 2, are shown in 
Fig. 2. The transmission probability T is plotted versus the gate potential Uo- Solid lines correspond to the RHF 
approximation. 
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FIG. 2. Plots of the transmission probability T versus the gate potential Uo for three different values of the barrier length, 
I — 8, I — 10, / = 12, and for three different values of the electron density in the "leads", no = 1, no = 1.5, and no = 2. Sohd 
lines show results of RHF calculations while dotted and dashed lines show spin up (Tf ) and spin down (T[ ) transmission 
probabilities calculated within the UHF method. Dotted and dashed lines in Fig. "a" (I = 8) are not distinguishable from solid 
ones. 

The RHF calculation for n w 2 has been performed earlier in Ref. 17 . All the plots presented in Fig. 2 clearly 
demonstrate structures of the conductance. An important point is that reduction of the electron density in leads 
and hence reduction of screening results in enhancement of the structure. Another feature is the evolution of the 



structure down for longer contacts. The results of UHF calculations are shown in the same Fig.l by clotted and 
dashed lines. The dotted line shows the transmission probability for the spin "up" channel and the dashed line 
shows the same for the spin "down" channel. Certainly the choice of "up" and "down" is arbitrary, one can swap 
the spins. The UHF method always gives two degenerate solutions. For I — 8 the UHF results are not presented 
because they are hardly distinguishable from that of the RHF method shown by solid lines. According to Eq. 
(2) the observable transmission coefficient is the average of Tj and Tj. The results of RHF and UHF methods 
are very close. To demonstrate the closeness we also present in Fig. 3 plots of electron densities n-\{x) and ni(x) 
for parameters (Uo = 6.1, n = 1.5) and (Uo — 6.4, n — 1.5) that correspond to two points on the shoulder in 
Fig. 2b. Solid lines represent the RHF density, nj(x) = n>i(x). The dotted line and the dashed line represent UHF 
densities n-\(x) and ni(x) correspondingly. In case (a) (Uo = 6.1, no = 1.5) the dotted and the dashed lines are 
not distinguishable from the solid one. In case (b) (Uo = 6.4, n = 1.5) they are distinguishable, but very close. 




FIG. 3. Electron densities n-\(x) and n±(x) for parameters (Uo = 6.1, no = 1.5) and (Uo = 6.4, no = 1.5) that correspond 
to two points on the shoulder in Fig. 2b. The solid lines show RHF density, n-\(x) = ni(x). The dotted line and the dashed 
line show UHF densities n-\(x) and n±(x) correspondingly. In the case (a), (Uo = 6.1, no = 1.5), the dotted and the dashed 
lines are not distinguishable from the solid one. 



According to Fig. 2 the RHF and UHF methods really disagree only at I = 12, n = 1, Uo > 3.7: relatively long 
contact, very low electron density in leads, and small conductance. This is the regime where the Kondo model is 
relevant. To understand what is going on in this situation we present in Fig. 4 plots of electron densities n^(x) and 
ni(x) at Uo = 4. As in the above figures, the solid line represents the RHF calculation, and dotted and dashed 
lines represent the UHF calculation. 
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FIG. 4. Electron densities for parameters I = 12, no = 1, and Uo = 4 that correspond to the structure in the left curve 
in Fig. 2c. The solid line shows the RHF density, n-\(x) = ni(x). The dotted line and the dashed line show UHF densities 
n-f (x) and (x) correspondingly. 

Clearly in this situation the RHF approximation is wrong. According to the UHF calculation the spin down electron 
density within the contact is practically zero and, on the other hand, there is one spin up electron localized in the 
contact. There is no doubt that in this case dynamics of the contact is Kondo-like. In this case even the UHF 
method does not give a correct transmission coefficient because the method does not take into account long-range 
Kondo dynamics. However, fortunately, the answer is well known; the transmission coefficient is peaked up to 
unity, see Ref. 20 . So, the correct plot of the transmission coefficient at I = 12 and n ~ 1 coincides with that 
presented in Fig. 2c for Uo < 3.7, and then there is a narrow peak up to T = 1 at Uo « 4. It is interesting to 
note that the transmission coefficients calculated within the Hartree-Fock approximation for no = 1 and for shorter 
contacts (Fig. 2a, b) have a qualitatively similar dependence: deep minimum and a narrow peak. This similarity 



clearly demonstrates how the several-electron bound state that can be assessed by the Hartree-Fock method (Fig. 2a, 
no = 1) evolves to the multi-electron Kondo bound state that cannot be assessed by this method (Fig. 2c, uq = 1). 

It is interesting that for longer contacts one can trap more than one electron in the contact. To illustrate this 
in Fig. 5 we show UHF electron densities n-f(x) and n^{x) for a contact of length / = 20 and density in the leads 
no = 0.56. In this case the two electron solution of the type shown in Fig. 5a is realized at the gate potential 
1.3 < Uo < 1.85, then at 1.85 < Uo < 2.15 the solution "jumps" to the single electron state shown in Fig. 5b. At 
the higher gate potential there arc no electrons in the contact. So adjusting the length of the contact, the density 
of electrons in leads, and the gate potential one can pin within the contact a single electron like it is shown in Fig. 4 
and Fig. 5b or even the two electron "molecule" shown in Fig. 5a. However, before getting to this very strongly 
correlated regime the transmission probability dips down to the value of few percent. At most the probability in 
the dip is 20% as it is shown in Fig. 2c (no = !)• There are no such dips in experimental data. Therefore it is 
unlikely that Kondo dynamics can be relevant to the effects observed in works 4-8 . On the other hand the plots 
shown in Figs.2a,b,c for tiq = 1.5 and uq = 2 look very similar to the experimental data. Structures on these plots 
are related to the few-electron spin dependent bound state. The closest physical analogy in this case is probably 
the Peierls spin-density instability. 
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FIG. 5. UHF electron densities n^(x) (dotted line) and ni(x) (dashed line) for a "very" long contact, I — 20, and for a 
"very" low electron density in leads, n — 0.56. 



In conclusion, within a one-dimensional model we have analyzed the conductance of a short quantum contact at 
zero temperature. Restricted (RHF) and unrestricted (UHF) Hartree-Fock methods have been used in the analysis. 
Both methods clearly demonstrate structures very similar to that observed in Refs. 4-8 . Agreement between RHF 
and UHF methods confirms the validity of both approximations. The conductance structure is related to the charge 
density wave developed in the contact. This is a spin dependent effect because without the exchange interaction 
the structure disappears, so this is a kind of spin-dependent bound state within the contact. Reduction of the 
electron density in the leads and hence reduction of the screening results in enhancement of the structure. The 
structure evolves down for longer contacts. 

Having the contact long enough, the density of electrons in the leads low enough, and adjusting the gate potential 
one can pin within the contact a single electron or even a two electron "molecule" . The single electron would imply 
Kondo-likc dynamics as has been suggested in Refs. 8 ' 18 . However, according to our calculations, before getting to 
this regime the transmission probability as a function of the gate voltage dips down to at least 20%. Such dip has 
never been observed experimentally. Therefore it is unlikely that Kondo dynamics can be relevant to the effects 
observed in works 4-8 . 
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